Classical spin systems and the quantum stabilizer formalism: general mappings and 

applications 
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We present general mappings between classical spin systems and quantum physics. More pre- 
cisely, we show how to express partition functions and correlation functions of arbitrary classical 
spin models as inner products between quantum stabilizer states and product states, thereby gen- 
eralizing mappings for some specific models established in our previous work [Phys. Rev. Lett. 98, 
117207 (2007)]. For Ising- and Potts-type models with and without external magnetic field, we show 
how the entanglement features of the corresponding stabilizer states are related to the interaction 
pattern of the classical model, while the choice of product states encodes the details of interaction. 
These mappings establish a link between the fields of classical statistical mechanics and quantum 
information theory, which we utilize to transfer techniques and methods developed in one field to 
gain insight into the other. For example, we use quantum information techniques to recover well 
known duality relations and local symmetries of classical models in a simple way, and provide new 
classical simulation methods to simulate certain types of classical spin models. We show that in this 
way all inhomogeneous models of g-dimensional spins with pairwise interaction pattern specified 
by a graph of bounded tree-width can be simulated efficiently. Finally, we show relations between 
classical spin models and measurement-based quantum computation. 
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I. INTRODUCTION 

Classical spin systems are widely studied in statistical 
physics [l|. They also play an important role in model- 
ing complex behavior also in other disciplines, such as 
economics and biology. In spite of their often simple 
definition, spin models show a highly non-trivial behav- 
ior, as is, e.g., apparent from their phase structure and 
criticality. Surprisingly, even the simple Ising model of 
interacting 2-state spins arranged on a 2D square lattice 
(with external magnetic field) is in general not solvable, 
and calculating, e.g., the ground state energy or the par- 
tition function is known to be a computationally hard 
problem 0|. 

In quantum information theory (QIT), on the other 
hand, (entanglement) properties of quantum systems 
are systematically studied, and possible applications re- 
garding, e.g., quantum computation are investigated. 
QIT has become a field of interdisciplinary interest, and 
concepts and methods developed in QIT have found ap- 
plications also in other branches of physics. In the con- 
text of QIT, methods to efficiently compute and simu- 
late certain quantum systems and their properties have 
been developed @J3 • In particular, so-called quantum 
"stabilizer states" (J (| 0] and "graph states" [8, §] have 
been introduced and studied in detail. Stabilizer states 
are used for certain types of quantum error correction [5j 
and measurement-based quantum computation [ltj , and 
can be described efficiently in terms of their stabilizing 
operators. This allows to determine many of their (en- 
tanglement) properties, and to efficiently simulate some 
processes classically. 

In this paper, we present general mappings between 



classical spin systems and quantum physics related to 
QIT. More precisely, we show how to express the par- 
tition function and correlation functions of an arbitrary 
classical spin system as a quantum mechanical ampli- 
tude (scalar product) between a stabilizer state \tp) en- 
coding the interaction pattern, and a certain product 
state ®j\otj) encoding the details of the interaction (i.e. 
the coupling strengths) and the temperature: 
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With such a mapping at hand, we can use methods and 
techniques established in one field to gain insight into 
the other, thereby providing a novel approach to these 
problems. We have initiated this approach in a recent 
publication [ll[, where such mappings have been estab- 
lished for Ising and Potts-type models. Here we gener- 
alize this approach, and discuss the mappings and their 
applications in more detail. 

We further note that connections between quantum 
information theory and statistical mechanics have re- 
cently been studied by several other researchers [H, [l3| . 



A. Mappings between classical spin systems and 
quantum physics 

In this section we briefly sketch the general form of 
the proposed mappings between classical and quantum 
systems. 

We consider classical q-state spin systems with an 
arbitrary pairwise interaction pattern, described by a 
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graph G with vertex set V (position of the classical 
spins) and edge set E (corresponding to interactions). 
Such systems are sometimes called "edge models" (i.e., 
the interactions take place on the edges). Each spin 
s may assume q different states: s £ {0, . . . ,q — 1}. 
We will consider models where the pairwise interactions 
h(s, s') between spins s and s' are of the following forms: 

(i) h(s, s') only depends on the difference (modulo q) 
of the two involved spins, h = h(\s — s'\ q ); 

(ii) h(s, s') is of the form (i), but with additional local 
magnetic fields; 

(iii) h(s,s') is completely arbitrary. 

We will also consider (iv) models with arbitrary fc-body 
interactions. 

The Ising- and Potts model without [with] magnetic 
field are of type (i) [(ii)] respectively, while so-called 
"vertex models" (i.e., the interactions take place on the 
vertices) are a special case of type (iv). 

In each of the cases (i)-(iv), we show how one can ex- 
press the partition function Zq as an overlap between a 
quantum stabilizer state and a complete product state, 
(Eq. |T|)). Depending on the different forms of the in- 
teraction (as in (i)-(iv)), these quantum states will be 
defined slightly differently. 

(i) For models without local fields, the correspond- 
ing quantum states consists of \E\ g-level quan- 
tum systems (one for each pairwise interaction 
term). We will denote the stabilizer state by \tpG)- 
The product state has the form \a) — ® eeE |a e ), 
where the coefficients of each \a e ) encode the 
strengths of the pairwise couplings, as well as the 
temperature of the system. 

(ii) For models with local magnetic fields, the corre- 
sponding quantum states consist of \V\ + \E\ q- 
level quantum systems (one for each pairwise in- 
teraction term and one for each local field), with 
stabilizer states denoted by \ipc} and a product 
state \a) = <g) e£E \a e ) <g) aeV \a a ). 

(iii) For models with general pairwise interaction (iii), 
we provide a mapping where the stabilizer state 
is a tensor product of \V\ entangled states, \4>) = 

(gWvG^ob'} 8 "")- Here > n * is the degree of 
vertex a, i.e. the number of neighbors in the graph, 
which also determines the number of associated q- 
level quantum systems. Correspondingly, we now 
consider states \a a b) of dimension q 2 for the over- 
lap, which are associated to one quantum particle 
belonging to vertex a and and one quantum parti- 
cle belonging to vertex b. A similar picture holds 
for models with arbitrary fc-body interactions (iv), 
where the product states have now dimension q k , 
and are associated with multiple vertices. 

We will investigate the entanglement properties of the 
states \"4>g) and \(po) and their relation to the underly- 
ing interaction pattern specified by the graph G, and 



provide a number of examples to illustrate this connec- 
tion. 

The mappings (ii)-(iv) can be extended, and will allow 
us to express also classical correlation functions in a 
quantum language. 

B. Applications of the mappings 

Based on these mappings, we will then illustrate some 
applications. Here we briefly sketch which applications 
can be obtained. 

(a) Using well established stabilizer methods @, @, 0, 
H], we show how one can recover the well known 
high-low temperature duality relations [H] for clas- 
sical spin models on arbitrary planar graphs. 

(b) Using the fact that stabilizer states are stabilized 
by certain tensor product operators, we derive lo- 
cal symmetry relations for classical models, i.e. we 
identify models with different coupling strengths 
that lead to the same partition function. 

(c) We show how one can use recently established re- 
sults in QIT to classically simulate certain classes 
of quantum systems efficiently 0, 0, EH and thus 
obtain novel simulation algorithms for classical 
spin system. More precisely, by describing sta- 
bilizer states in terms of an optimal tree tensor 
network [j| of dimension d, one can compute the 
overlap with product states with an effort that is 
polynomial in d. This leads to an efficient algo- 
rithm to classically simulate arbitrary (inhomoge- 
neous) classical g-state models on graphs with a 
bounded (or logarithmically growing) tree width. 
We also extend these results to models with fc- 
body interaction. 

(d) Finally, we discuss links between classical spin 
models and measurement based quantum compu- 
tation. This allows us to relate the computa- 
tional complexity of computing partition functions 
of classical spin models with the quantum compu- 
tational power of the associated graph states. 

We also note that (d) has recently been used in 
Ref. [TEI to show a " completeness" property of the 
2D Ising model. That is, invoking the connection 
to measurement-based quantum computation, it was 
shown that the partition function of any model with 
pairwise interaction in arbitrary dimension can be ex- 
pressed as a special instance of the partition function 
of a 2D Ising model on an (enlarged) 2D square lattice 
(with complex coupling strengths). 

C. Guideline through the paper 

The paper is organized as follows. We start in Sec.HTl 
by briefly reviewing classical spin models, and collect 
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some relevant results on stabilizer and graph states in 
Sec. IIIII We then introduce different mappings between 
classical spin systems and quantum mechanical ampli- 
tudes, and discuss the properties of the involved quan- 
tum states in Sees. HVl and fVl We illustrate a number 
of applications of these mappings in Sec. IVI[ and sum- 
marize and conclude in Sec. IVIII 



II. BACKGROUND ON SPIN MODELS 

In this section we describe the classical models that 
we want to consider. Since the various approaches to be 
described later are related and can be viewed as deriva- 
tions from an original scheme, we will focus on the orig- 
inal approach first. 

The typical model to be considered by the original 
approach is the thermal state of a classical spin model 
described by a Hamiltonian function with two-body in- 
teraction, and this model will serve as an introductory 
guide to the general idea. These systems have the virtue 
that they admit a description by means of a graph : 
the spins of the system correspond to the vertices and 
the two-body interaction pattern between the spins is 
given by the edge set. 

We will describe a mapping of such an interaction 
graph to a stabilizer state of a quantum system. Per- 
forming an overlap of this quantum stabilizer state with 
another quantum product state, encoding the temper- 
ature and individual interaction strengths, then yields 
the properties of the thermal state of the classical sys- 
tem. We want to emphasize that this evaluation is not 
approximate but exact. Later on, extensions of this for- 
malism will be given as well, going beyond this partic- 
ular kind of graphical description and at the same time 
going beyond the limitation to two-body interactions. 

It is important to keep in mind that the interaction 
pattern and the interaction strengths are encoded at dif- 
ferent places: the graph encodes the interaction pattern, 
not the strengths, hence an edge connecting two vertices 
simply denotes the fact that there is an interaction tak- 
ing place. The strength and nature of this interaction 
is not encoded in the graph, but in a product state to 
be specified later. This encoding admits the strengths 
of all edge terms and all vertex terms to be chosen in- 
dividually, hence the interaction strength may vary for 
different pairs of spins and also the local field may vary. 

More precisely, let, for now, H be a Hamiltonian func- 
tion with two-body-interaction between classical spins s 
that can assume q possible values s e {0, —,q — 1}. In 
the graphical description of this Hamiltonian function, 
we let G = (V,E) denote the graph associated with H, 
where the sets V and E contain the vertices and the 
edges of the graph respectively. In this picture, any ver- 
tex v € V corresponds to a classical spin site s v and any 
edge e € E between adjacent vertices v%, v 2 of the graph 
corresponds to an interaction term between the respec- 
tive spins s vi and s V2 . Additionally, we allow each spin 
s v to contribute a local term to the Hamiltonian func- 



tion, i.e. a term that that depends on the state of the 
site s v alone, although this is not reflected in the graph. 
We might think of the energy of the spin in a local field. 
We choose the graph to be a directed one, denoting the 
orientation by a. The exact choice of the directions can 
be arbitrary but has to be fixed. This way, the two ad- 
jacent vertices of an edge e £ E can be distinguished as 
"head" v+ and "tail" v~ of the edge, respectively. 

We will derive several different mappings for Hamil- 
tonian functions described by these graphs. The first 
mapping admits descriptions of systems with classical 
Hamiltonian functions of the form 

= £Mk. + -*«.-!,)' (2) 

with h e being an energy term that depends on the rel- 
ative state of two interacting spins s v + and s - modulo 
q. In the second mapping we extend the quantum de- 
scription to be able to include also external fields 

H {{Si}) = £ fee (| V ~ s v - \ q ) + J2 & "M> ( 3 ) 

eeE v£V 

where b v is an energy term contributed by a local ex- 
ternal field, acting on the spin s v . To go beyond the 
limitation to interaction Hamiltonian functions that de- 
pend on the relative state of the spins only, we finally 
provide further approaches to treat Hamiltonian func- 
tions of the form 

H({si}) = £ %j')(si,Sj) 

as well as arbitrary Hamiltonian functions with ro-body 
terms. 

The degrees of freedom in the definitions of these 
Hamiltonian functions give rise to a large set of pos- 
sible classical spin systems to be described - even if 
we restricted ourselves to the sets of Hamiltonian func- 
tions specified in Eqns. (J2j) and ((SJ. Among those are 
the Ising model, the Potts model and the clock model 
on arbitrary lattices, all equipped with (local) magnetic 
fields, and generalizations thereof 

A. Ising model 

The Ising model describes a set of classical spins (or 
simply dipoles) that can point either up or down and are 
placed on a graph. All next neighbors have the same dis- 
tance (hence the interaction strength is uniformly given 
by the real number J) and long-range forces are ne- 
glected. Moreover, there is a global external field whose 
strength is given by the real number B, which puts an 
energetic bias on the possible configurations. Thus the 
classical Ising model is described by the Hamiltonian 
function 

H lsins ({*}) = JJ2 h - s j\ 2 - \)> (4) 
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where the S{ £ {0, 1} and denotes that i and j are 
adjacent spins on the graph. We note that it can be 
rewritten as 



Hr^ s (W) = -J' E 



2J'5> 



where cr, G {+1,-1}. This is the more familiar form 
and can be obtained from Eq. H| by a rescaling of 
parameters and an addition of a constant. Although 
this model is highly idealized, it features (in appropri- 
ate dimensions) many properties of realistic solids, such 
as phase transitions, spontaneous symmetry breaking 
etc. As will be shown, our treatment allows — without 
a change of computational effort — the generalization to 
spin-glass Hamiltonian functions, where the factor J is 
actually dependent on the specific pairs of spins that 
interact: J — > Jy. 



B. Potts and clock models 

A generalization of the Ising model is given by the 
Potts- and the clock model. Whereas the individual 
spins in the Ising model can take only one of two values 
and hence for neighbors there are only the alternatives 
of being parallel or anti-parallel, it might be desirable 
to allow the individual dipoles to assume more posi- 
tions and hence to obtain more relative configurations of 
neighbors that can be discriminated energetically. Ac- 
cordingly we choose spin states s< <E {0, q — 1} and a 
Hamiltonian function 



H({si}) 



E J & 



(5) 



where 9^ is a function that discriminates the relative 
states of neighboring spins. We can interpret it for 
instance as the angle between adjacent dipoles, pro- 
vided that they can only rotate in a fixed plane, e.g., 
©ij = 9j — @j with discretised positions &i = 2-rrSi/q. 
The function J, which characterizes the Hamiltonian 
function, maps the relative angle (i.e., relative state) of 
adjacent spins to an energy value: The Potts model is 
defined by 

•WO,,) :=-eS(e. l3 ) 
with e£l and the clock model by 

^dock(0y) := ~ecos(Q lj ). 



central quantities of interest that we want to obtain are 
the partition function 



Z{f3) 



5> 



as well as the n-point correlation functions, whose defi- 
nition can be found, e.g., in Ref. [l7| 
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= Z cos (©ij cos (6j 3 ) ... cos (©,„) e 



The partition function encodes the macroscopic prop- 
erties of a thermal ensemble. The parameters that en- 
ter depend on the kind of ensemble we look at, e.g., 
the canonical (temperature), grand canonical (temper- 
ature and chemical potential) and others. In the present 
framework we will deal with the canonical ensemble, be- 
cause the number of spin sites is fixed, but energy can 
be drawn from an external bath. 

Let us briefly illustrate the importance of the par- 
tition function. The partition function of a canonical 
ensemble is 



E« 



-0Et 



where the index i is the index for the states with energy 
Ei that the system can take and (3 — (fcsT) with the 
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Boltzmann constant fcs- Moreover, pi = Z~ e 
the probability to find the system in the state with en- 
ergy Ei . Several relevant quantities can now be derived 
from Z: We can extract the expectation value of the 
energy 



(E), 



Z 



-pEi 



d log Z 

d(3 ' 



the variance of the expected energy 

d 2 \ogZ 



(5EY 







d(3 2 



as well as the free energy 



F={E) fj -TS = -f3- 1 \ogZ, 

where the entropy is S = — fcs J2i Pi l°gPij an d more. 
We refer the reader to standard text books on this topic. 



III. STABILIZER STATES AND GRAPH 
STATES 



C. Partition function 

The focus of this paper will be on the thermal equi- 
librium of these classical systems. More precisely, the 



In this section, we give the definition and some prop- 
erties of stabilizer states 0, 0, 0] and graph states 0, Q . 
We will first consider spin- 1/2 quantum systems, then 
proceed to higher dimensional systems. 
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A. Graph states 

Here we will briefly familiarize the reader with the 
graph states. In the present context, a graph G = (V, E) 
is identified with a quantum system. Each vertex a rep- 
resents a quantum spin, and the adjacent vertices (con- 
nected with a by edges in the graph) form the neigh- 
borhood iV a of a. This way, the graph defines a set of 
operators 



K a := J] 



b£N a 



where the sigma-matrices are defined as usual 



cr 



I 

1 

-% 

1 



, Ox 



1 

1 

1 

-1 



(6) 



and the notation 0' a ' of an operator O means the tensor 
product of the operator O, acting on the subspace of site 
a, and 1 everywhere else. A graph state \G) associated 
with to the graph G, and hence with the set {K a }, is 
the unique non-trivial fixed point of the operators K a , 

VK a : K a \G) = \G) . 

Graph states are a subset of the stabilizer states, 
which play an important role in the context of one-way 
quantum computing. Conversely, every stabilizer state 
can be written, up to a local rotation, as a graph state. 



B. Stabilizer states 

We will now turn our attention to the slightly more 
general set of stabilizer states. The concept of defining a 
state as a simultaneous fixed point of a set of operators 
can be used in a slightly more general way than in the 
case of graph states, where the operators K a take a 
very special form. To construct more general sets of 
operators we consider the sigma-matrices, see formula 
([6]), and the group they generate 

Qi = {±a ,±icro,±cr x ,±icr x ,±ay,±icry,±a z ,±ia z } . 

Tensor products of Q\ with itself form the Pauli groups 
Qn '■— Qf n - It is known that any Abelian subgroup 
S C Q n of a Pauli group with |«S| = 2" that does not 
contain — l n has a unique fixed point in the Hilbert 
space H that it acts upon. We then call S the stabilizer 
of and \ip) a stabilizer state. It should be noted 
that each stabilizer can be identified with its generator, 
i.e., a set of operators that generate it. Generators are 
not unique sets, but share the necessary requirement to 
contain n independent operators. 

For our purposes, the prefactor (±1, ±i) of an element 
of a Pauli group will not be important. Moreover, there 



is a mapping between the Pauli group Q n / ~ {Qn modulo 
prefactors) and the group F 2 .™, which will be used later. 
Since <r y = icr x a z and <7° = I2 for all sigma-matrices, 
we can encode the generators of Q\j ~ as follows 



cr - 




-(00) 


Ox r 




-(10) 


Oy - 




-(11) 


e z r 


- ->l - 


-(01) 



where ~ denotes equality modulo prefactor. Tensor 
products of these operators and hence elements of the 
groups Q n / ~ will be encoded by the mapping 

n 
i=l 

The generalization to g-dimensional quantum systems 
with H = (C 2 )® 9 is straightforward. We replace a x and 
a z by the operators X and Z respectively, where 



X \j) =\j + l mod q) , Z \j) = e 



2mj/q 



q — 2 being a special case that gives us back a x and a z . 
The higher-dimensional groups Q^l ^ are thus generated 
by tensor products of X a Z b where a, b = 0, q— 1. The 
mapping is generalized to the group homomorphism 



(0«/~,-)3®**Z C< 
i=l 

~ (6,...,^,Cl-,Cn) G (Ff ,+). 

The number of elements in a stabilizer that stabilizes 
one single stabilizer state is q n , the number of elements 
of its generator is n. 

Related to this construction is a theorem that we will 
use later. Note that we do not neglect the phase this 
time. 

Lemma 1. Any two operators ^"^X^Z^ and 
<2)™_i X^'Z^ commute if and only if £' • £ — £ • £' = 
modulo q. 

Proof. The computation for the single spin site yields 

= X £iZ^Xt*Z^e 2 ™^ i -^y«. 
Hence for all sites together we obtain a phase factor 

e 27ri(£'<-£-e')/9. □ 

It is noteworthy that for q — 2 each stabilizer state is 
related to a graph state by some local unitary transfor- 
mations. This means that the two sets do not differ as 
far as their non-local properties are concerned. 

The stabilizer states are interesting to us, because — as 
will be shown — the interaction patterns of the Hamil- 
tonian functions of the classical spin systems that we 
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look at correspond to such states. Moreover, stabilizer 
states are well investigated and elaborate techniques are 
known for their manipulation [5], allowing us to investi- 
gate relationships between different (interaction) graphs 
and hence different Hamiltonian functions. 



IV. ENCODING CLASSICAL SPIN SYSTEMS 
IN QUANTUM LANGUAGE 

In this section we will investigate in detail the cor- 
respondence of the classical and the quantum systems 
that were presented in the preceding sections. 



A. The basic principle 

The basic approach, which was introduced in 
Ref. [11], is sufficient to describe systems with clas- 
sical Hamiltonian functions of the form H ({sj}) = 



The idea is to map the graph 



G, describing the interaction pattern into a stabilizer 
state, together with a supplementary product state that 
encodes the interaction strengths as well as the temper- 
ature. 

Let the classical spin system be defined by the (arbi- 
trarily oriented) interaction graph G a — (V, E) over | V| 
classical spins of dimension q, where a denotes the ori- 
entation. Let in the following M = \V\ and N = \E\. 
Now consider the incidence matrix B a of the interac- 
tion graph G a . This matrix has one row for each vertex 
and one column for each edge. The entries are either 
or ±1, where B° e — —1 if the index pair (v,e) corre- 
sponds to the tail vertex v of edge e, e = +1 for the 
head vertex v of edge e and B° e — otherwise. Consis- 
tent with our notation, we do not consider graphs with 
edges that connect one point with itself. The rows of 
B a span the Z 9 -vector space Cg (q), which is a linear 
subspace of . The vectors (B a ) T s £ Cg {q), where 
T denotes transposition, correspond to the vectors that 
encode spin configurations (s v ) veV , as the linear map- 
ping 



v£V 



shows. 

Lemma 2. The kernel of the linear mapping (B a ) T 
has q K elements, where k is the number of connected 
sub-graphs of G (without isolated points). 

Proof. We re-arrange the rows of the matrix of (B a ) T 
so that the connected sub-graphs Gi — (Ei,Vi) are de- 
scribed by blocks Bi, i.e., 



( B x 

B 2 

B 3 

\ i ; : 



Within each connected sub-graph Gi, there is at least 
one path from each vertex v to each other vertex v': 
(v, Vq, v\, v'), each edge (v n ,v n +i) in this path be- 
ing represented by one row in the corresponding matrix 
Bi. Since a vector s to be in the kernel of Bi implies 
\s Vn — s v n+1 \ q — for each edge (v n ,v n+ i), we deduce 
immediately that \s v — s v >\ q — for any two vertices v 
and v' inVi. Hence if s is in the kernel of Bi , all spins in 
{ s v n } v e y. take the same value. So there are q different 
vectors in the kernel of each matrix B^ , of which there 
are k. □ 

We are now ready to define an non-normalized sta- 
bilizer state encoding G° '. We obtain it by first in- 
terpreting each vector c = (c±, C2, cat) € Cg (q) as 
a product state of a multipartite quantum spin sys- 
tem with spin-dimensionality q according to the formula 
\c) := |ci) ® |ca) <8> ... <8> Icy) and by a subsequent sum- 
mation of all these states [23] 



q 



cec G (?) 



(7) 



where the second equality and the factor q K follow im- 
mediately from lemma ([2]). For an illustrative example 
see Fig. [TJ 

Lemma 3. The state \4>g) is a stabilizer state. Its sta- 
bilizer consists of the q N operators 

X{v)Z{u):=§§X*°Z u °, (8) 

where v £ Cg (q) and u £ Cg (<?) ■ 

Proof. From Lemma (JXJ we derive immediately that, by 
using the given construction rule for the operators, we 
obtain a commuting set. Considering the equation 

and hence 

X (v) Z (u) X («') Z (u') =X(v + v')Z (u + u') e 2 ™ u - v '/«, 

where u ■ v' = J2i u i v [ = f° r each admissible choice of 
these vectors, we also see that these operators form a 
group. Furthermore, these operators actually stabilize 
the (nontrivial) state \iPg), since for all v £ Cq (q) and 
for all u £ Cg {q) 

X(v)Z{u)\i> G )=X(v)Z{u)q* J2 l c > 

cGC G (g) 
= q K J2 X(v)Z{u)\c) 



c£C G (q) 
c'eCa(q) 



e 2mu-c/ q ,_ 
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edge of 

interaction graph G classical 
, M ^ site 




quantum 
site 



FIG. 1: The basic construction principle. This figure shows an example of an encoding of a classical interaction pattern into 
a stabilizer state. Thin graph: the classical interaction graph G; thick graph: the derived graph relating quantum sites in a 
stabilizer state. The classical spin sites correspond to vertices in a graph G. The interacting pairs of sites are mapped to a 
quantum site, one for each edge ("edge qudits"). The quantum sites form, by construction, a stabilizer state. 



From this, we can moreover deduce that —1 is not ele- 
ment of this set of operators. 

For this set to be a stabilizer of a single state of our 
Hilbert space, the number of elements in this set has 
to be q N . This part of the proof is given in Appendix 
EJ □ 



1. Thermal quantities 

Now we are able to formulate the central theorem of 
this section. 

Theorem 4. The partition function Zq (q, {h e }) of a 
classical spin system defined on the graph G = (V, E) by 

„+ - 



the Hamiltonian function H ({sj}) = ^2 eeE h e (\s v -t 
s v - | ) can be written as the overlap of a stabilizer state 
and a product state 

Z G (q,{h e }) = ((g) (OeDl^o) 
eS-E 

of a quantum mechanical spin-system, where 



a,. 



9-1 

£« 

i=o 



-0KU) 



li) 



(9) 



Proof. The state \ipa) ls a stabilizer state according to 
lemma |(3]), and we compute, with an arbitrarily chosen 
orientation a of the graph G, 



(®<«e|)llM = £((g)<ae|) (B ff fs 

(9) 



-0hMs + -s _ ) 



which concludes the proof. 



Let us give a brief interpretation of the method used 
to encode the partition function. We observe that to cal- 
culate partition functions of systems with Hamiltonian 
functions of the form 



eeE 



(zero external field) it is already sufficient to know the 
relative state of spins whose corresponding vertices are 
connected by an edge. Accordingly, we map each vec- 
tor (s v ) veV of spin configurations to the corresponding 



one (B a ) 



S„+ - 8„- 



of differences along 



"7/eS-E 

edges using the incidence matrix B a . These vectors are 
automatically consistent with spin-configurations, and 
moreover, there can be no more of them than we have 
already given. 

As shown, an interaction pattern is encoded into a 
graph and this graph is encoded into a stabilizer state. 
Furthermore, the corresponding interaction strengths 
(as well as a temperature) are encoded into a product 
state. This way we encode all the information about 
the partition function of a thermal state into two states 
with comparatively simple structure. 

Example 5. Here we consider examples of states \a), 
which encode the interaction strengths of the Hamilto- 
nian function. For the models we consider, these are 
product states \a) = <S) e e_B \ a e), which are derived im- 
mediately from the respective Hamiltonian functions 
given in section |TT|) . 

1. For the g-state Potts model the state \a) is de- 
rived from the Hamiltonian function ©, with the 



function J given by </ Pot t B (Si 



-e5(Qij). This 



Hamiltonian function is characterized by two- 
body interactions, whose strengths are encoded 
into states \a e ) which take the form 



= e" B 



□ 



q-l 
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2. For the g-state clock model the state \a) is derived 
from the Hamiltonian function |(5]) , with the func- 
tion J given by J c i ock (Oy) := — ecos(Oy). The 
individual two-body interaction strengths are thus 
encoded into states 

9-1 

l«e)=|«) clock =E e/3£COS( ^ /9) tf>" 

3. As a special case, for q = 2 we obtain, in an anal- 
ogous fashion, the states \a) and |a e ) for the Ising 
model 

\a e ) = \a) lsing = \0)+e-P J \l). 

In the following part we look at examples of states 
\iPg), which encode the interaction patterns of associ- 
ated Hamiltonian functions, thereby investigating spe- 
cial cases of graphs and their corresponding stabilizer 
states. 

1. Tree Graphs. Here we consider models whose 
interaction patterns are characterized by tree 
graphs, i.e., graphs containing no loops. The 
statement that "n columns {cj} i=1 of the inci- 
dence matrix B a of a graph G are linearly depen- 
dent" means that there is a non-trivial linear com- 
bination such that Y17=i ^ iC i = 0- Hence there is 
at least one vector that equals the negative sum 
of the remaining ones, say, c\ = —A] -1 ^T=2 ^ c «- 
Since the columns describe the start and end 
points of the edges, this means that the graph 
contains a loop. In turn, loop-less graphs (= tree 
graphs) have an incidence matrix with N = \E\ 
linearly independent columns and hence N lin- 
early independent rows. This means that the rows 
span the entire space (= Cg (?)) and hence 

Eli) • (io) 

j=o J 

In conclusion, we observe that the states derived 
from tree-graphs are product states. 

2. A cycle. Here we consider models whose inter- 
action patters are cycles, i.e. a closed loop. If the 
graph is a closed chain, the incidence matrix looks 
(besides reordering of the edges) like this 

/-ll \ 

-1 1 ••• 

B <y = 0-1 

: '•• 1 

V 1 ■■■ -1 / 

We see that a vector v that is perpendicular to 
all rows has the property = Vj for all i and 



j, and hence Cq {q) ± = span |(1, 1, 1, 1, 1) T j- 
We hence choose 

{Z* N ,XW (X- 1 )^ \n = l,...,N-l\ 

as generating set of the stabilizer. We can verify 
that the state |^ G ) = (\jx)® N ^J , where \j x ) 

is an eigenstate of the A-operator, is an eigenstate 
of the generator of the stabilizer and hence the sta- 
bilizer itself. This state is invariant under reorder- 
ing of the edges and hence the proof is independent 
of the choice of B a that was chosen in the begin- 
ning. Thus, the states derived from graphs that 
are closed chains are (generalized) Greenberger- 
Horne-Zeilinger states (GHZ states.) In particu- 
lar, for q = 2 one obtains the state 1+)®^ + !-)®^ 
(where |+) and |— ) are the eigenstates of the Pauli 
matrix a x ). 

3. The Kitaev model. The Kitaev model of topo- 
logical^ protected quantum states is defined as 
follows. On each edge of a toric lattice with 
checkerboard structure we place one qubit, the 
edge qubit. The toric code state (actually a sub- 
space) is the common eigenstate of a set of opera- 
tors that are constructed using the neighborhood 
relations of the toric lattice. More precisely, for 
each but one of the smallest possible loops Li (the 
plaquettes) in the lattice, we define one operator 

b 1 ■= Yl z{a ' b) - 

(a,b)£Li 

We leave out one because it would not be inde- 
pendent from the others. Similarly, each vertex a 
(there is no qubit in the vertices) has a neighbor- 
hood N a of adjacent edges, forming a star. On the 
qubits of each but one of these stars we define the 
operators 

A a := J] X^ b \ 

One has to be left out because it is not indepen- 
dent of the others, as in the case of the plaque- 
ttes. All these operators mutually commute, be- 
cause in each loop a vertex has either zero or two 
nearest neighbors. Hence, these operators gener- 
ate a stabilizer, whose fixed point is the toric code 
state. We notice that this stabilizer consists of 
2 2N ~ 2 independent operators defined on a 2iV-site 
quantum system and hence the stabilized object 
is not a single state but a subspace of dimension 
4. We remark that the connection between the 
2D Ising model and planar (toric) code states was 
first proven and utilized in Ref. [13| . 
In view of the huge variety of classical spin mod- 
els and their interaction graphs, we want to point 
out that this state can be defined more abstractly 



9 



and more closely related to the -B^-matrix con- 
struction used in the other examples, as shown in 
the following: We assume q = 2 and consider an 
arbitrary graph G a = (V, E) with the essential 
property to contain N — 1 = |V| — 1 independent 
loops {Li} i J[ 1 . The loops now naturally define a 
specific neighborhood N a of each vertex a, namely 
the union of the sets of nearest neighbors of a in 
each loop N a = \J L . {b e V; (a, b) G L,}. With 
the loops and neighborhoods specified, we define 
as above the operators 



(a,b). 



B, 



beN a 



n zia ' h) - 

{a,b)£Li 



All of these operators mutually commute, because 
in each loop a vertex has either zero or two near- 
est neighbors. There are N — 1 independent oper- 
ators A a and N — 1 independent operators Bi, 
which can be seen as follows. Considering the 
operators Bi, the statement Yl ieI Bi = 1 with a 
set of loops I implies that I contains dependent 
loops, which is impossible for |/| < N. Similarly, 
considering the operators A a , for each set of ver- 
tices V the identity JlaGV" A a = 1 means that 
the sets N a — {(a, 6), b £ N a } (when consid- 
ered together) contain each edge twice. This is 
impossible if \V \ < N by construction of the in- 
teraction graph. On the other hand Jlaey A a = 1 
and \\ i Bi = t by similar arguments. As a special 
instance we recover the example given above if G 
is the periodic two-dimensional lattice, where the 
common fixed point of the operators A a and Bi 
defines the toric code state fl8j |. as introduced in 
the context of topological quantum computation. 

We will use yet another generalization of this con- 
struction procedure in the subsequent section, in order 
to access correlation functions and partition functions 
of systems with external magnetic fields. 



B. External fields and correlation functions 

The encoding scheme discussed in the last section is 
neither suited to evaluate correlation functions nor par- 
tition functions of systems with external fields, like those 
described by 

H ({Si}) = h e(\ S v+ ~ V L) + E b -( S -)- 



eeE 



To overcome this limitation we (have to) use a different 
encoding scheme. Instead of the state \ipa) we will now 
use the state 



E 

sG 



| S )|(S-) T s) 



N 



to encode the interaction pattern, where B a is again the 
incidence matrix of the interaction graph G. 



Lemma 6. The state \<pg) is a stabilizer state. Its sta- 
bilizer is generated by the N operators 



K a = n 

e:3beVs.th.(a,b)=eeE 



for every a 6 V and for every e = (a, b) € E, where 
a is either +1 or —1, depending on the orientation of 
the edge (a := B° a = —B a eh ). In our notation, an 
upper index in brackets denotes the audit acted on by 
the operator. 

Proof. We have to show that these operators have 
\<Pg) as a fixed point and, since the number of operators 
defined this way is TV = |V| + \E\ and hence equals the 
number of qudits in the quantum system in state \ fa), 
we have to show that they are independent. Under this 
condition they generate the stabilizer of a single state. 
To see the stabilizing property of the operators K a , we 
compute 

X {a) \{ (X^Y \s)\(B a ) T s) 

e:3beVs.th.(a,b)=eeE 

= \ S i)\(B°) T S >) 

with s' = (so, s a + 1 mod q, s \v\)) because 
X^ |s> = |s') and 

n (x^y^fs) 

b:(a,b)=e£E 



(B a fs) modq)= {B a fs'), 



where c a is the ath column of B a . Likewise, K e \ 
\<Pg), because 

(z^y° (z {b) y\s) 

= ( z ^y B - ( z ^y Blb is) 

= exp{-2™ (5 e > a + Bl b s b ) /q} |(S CT ) T s), 



and 

Z^\{B a ) T s) 

= exp{27ri (Bl a s a + B° b s b ) /q} \{B°) T s), 

so the phases cancel. 

The set of 2n operators that were just defined are 
mapped, by the isomorphism F^ n , to a set of 2n vectors 
that can be arranged in the following matrix 

( Vi T o \ 

{B a ) 
— B a 

V o t m j 



10 



This matrix has full rank, considering the Is. Hence 
the operators generate the full stabilizer. □ 
In the case q = 2 we recover a true graph state by 
an application of Hadamard transformation on the edge 
qubits. For an illustrative example, see Fig. [2j As be- 
fore, the classical spin sites correspond to vertices in the 
interaction graph G of the classical model. The inter- 
acting pairs of sites are then mapped to a quantum site, 
one for each edge (the edge qudits). What is different 
from the original scheme is that the individual classi- 
cal spin sites -acted on by local fields - are mapped to 
quantum sites as well, one for each vertex (the vertex 
qudits). The resulting graph is called a decorated graph. 
The resulting many body quantum states are again, by 
construction, stabilizer states. 



1. Thermal quantities 

We now come to the central result of this section. 
By means of the state \<Pg) an d appropriately chosen 
product states, we can compute the partition function 
of systems with local external fields as well as n-point 
functions. 

Theorem 7. The partition function Zq ({h ei b v } , 0) of 
a classical spin system at inverse temperature /3, defined 
on the graph G = (V, E) by the Hamiltonian function 
H ({si}) =EeGE /l e(|s w + - s v - \ q ) + J2vev b v ( s ") , can 
be written as the overlap of a stabilizer state and a prod- 
uct state 



Z G ({h e ,b v },{3) = 



((g)«l(g)(«e|)b G ), 

veV eeE 



whe 



\J) 



a, 



9-1 

E e 

3=0 
3=0 



Proof. The state \<Pg) ls a stabilizer state according to 
lemma |(6|), and we compute, with an arbitrarily chosen 
orientation a of the graph G, 



Likewise, we can write down a theorem for the n-point 
correlation functions. 

Theorem 8. The n-point correlation func- 
tions (si 1 ,Si 2 ,...,Si n )g of a classical spin sys- 
tem at inverse temperature (3, defined on the 
graph G a = (V, E) by the Hamiltonian function 
H{{si})=Y,e&E h e{\s v + ~ s v - \ q ) + T,vev h v (*«) , can 
be written as an overlap of a stabilizer state and a 
product state (up to a factor of Z , which is the partition 
function of the classical spin system). More precisely, 

(s il ,s i2 ,...,Si n >^ = Z- 1 (§Q(a' v (ii,...,i„)\§Q(a. e \)\\<pG) 

veV eeE 



where 



9-1 

E< 

3=0 
9-1 



x' v (ix,...,i n )) = ^cos(27rj/g) 



3=0 



and m v is the number of occurrences of v in the n-tuple 

(ix, ....in). 

Proof. The state \<p G ) is a stabilizer state according to 
lemma ([6]), and we compute 



((g) (a' v (ix,..., i»)|®(ae|)|VG> 
vev eeE 

= ((g)K(ii,...,i„)|(g){a e |) I s ) i^f* 

vev eeE 



£ n«»(2W9re-'" ,( " ) II 

sez«tiev eeE 

= cos(e n )cos(e i2 )...cos(e 4 Je^^« s '» 



where = 2irSi/q. We compare this to the definition 
of the n-point correlation function given in section UU 
This concludes the theorem. □ 



EXTENDING THE FORMALISM 



((g)«|(g)K|)| w ) 

vev eeE 

= ((g)K|(g)<a e |) J2 I s ) 

vev eeE seW 



(B°f a 



e n< 

sez™ vev 



-f3K(s v ) 



—j3h e (| S + — S 

= E^ 

sGZ" 



-em**}) 



□ 



A. The most general framework 

So far we have used product states of single edge-qudit 
sites, namely the states \a) — ® eeE |a e ), in the overlap 
with the states \<Pg) an d IV'g), to calculate partition 
functions and correlation functions. Allowing for tensor 
products of entangled states, \a) = ® eCE |«e), where 
the e are subsets of E with few elements, extends the 
set of possible encodings of classical spin systems. This 
is the content of this section. 

One shortcoming of the encoding of the interaction 
graph into the states \ and \^g) is the inability of the 
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edge of 
interaction graph G 



classical 
site 




quantum 
site 
"vertex qubit" 



|s,> 



quantum 

site 
edge qubit" 

Mi) 



FIG. 2: The extended construction principle. This figure shows an example of the extended encoding of a classical interaction 
pattern into a stabilizer state. Thin graph on the left: the classical interaction graph G; thick graph on the right: the derived 
graph relating quantum sites in a stabilizer state. The classical spin sites correspond to vertices in a graph G. The interacting 
pairs of sites are mapped to a quantum site, one for each edge (edge qudits). The individual classical spin sites on which the 
local fields act are then, too, mapped to quantum sites, one for each vertex (vertex qudits) - this is different from the original 
scheme. The resulting graph is called a decorated graph. The quantum sites are, by construction, in a stabilizer state. 



interaction to distinguish between classical spin states 
that have the same relative state \sa — s,-| = Is' — s'A 

1 J l q I L J I q 

but have different values Si ^ s' { , Sj ^ s'j. This in- 
ability stems from the fact that an attempt to encode 
pairs of neighboring states (si,Sj) into one edge qudit 
(si,Sj) i — ► \eij) via the _B-matrix formalism does not 
lead to a stabilizer state and hence fails, if |ey) takes 
more states than each of the sites Si or Sj . One way out 
of this dilemma is to encode the pairs of neighboring 
spin sites in the graph of the classical model into more 
than one qubit, while extending the overlap state \a) 
to states beyond product states. Although these states 
are not product states anymore, we can still interpret 
them as product states of composite particles, extend- 
ing over few sites as we restrict ourselves to subsets e 
of E with few elements. The entangled states moreover 
include neighboring sites only, which adds to the picture 
of composite sites (quasi-local states). 

We describe this generalization now and investigate 
the relationship to the more specialized cases. Under 
certain assumptions concerning the classical Hamilto- 
nian function, a formal mapping from the most general 
case to the more specialized ones is possible. Taking this 
step, i.e. performing this formal transformation, gives 
us a mathematical picture which is often much more 
enlightening than the original one. 



B. Encoding m-body interactions, each site 
appearing in maximally n terms of the Hamiltonian 
function. 

The most general case to consider is the one where 
we 



• allow each classical spin to appear in as many as 
n terms of the Hamiltonian function 

• allow each site to interact with m — 1 others 
(Hamiltonian function with m-body terms) 

• allow all configurations of the m interacting spins 
in each term to be differentiated energetically. 

Note however, that a simulation of thermal states of 
these systems on a classical computer scale unfavorably 
in m and n, as we will see in section IVI CI 

The first point in the list is addressed in the following 
way. Since each site is allowed to take part in n in- 
teractions, we need n instances of it in the stabilizer 
state. Of course, all instances of the local quantum 
systems have to be in the same state when measured. 
Hence we map each site to an n-body GHZ state: 
ei i— ► \ s i,i s i-2--- s i,n) ■ To address the latter two points 
of this list, we consider the following. To create a quan- 
tum state |7g) that enables us to differentiate energet- 
ically between all possible spin configurations of an m- 
body interaction, we map each site e, taking part in the 
interaction to a single quantum spin state i— > |sj) e .. 
The corresponding state \a), which is used for the con- 
traction that yields the partition function and which in 
the preceding sections used to be a product state, conse- 
quently has to be an entangled state in this picture. On 
the sites {e^}™ 1 taking part in one m-body interaction, 
the state \a) takes the form 

\a)= J2 e-^' 1 '"--^ \s lS2 ...s m ) . 

(si,S2, — ,S m ) 

Note however, that a simulation of thermal states of 
these systems on a classical computer scales unfavor- 
ably in m and n, as we will see in section IVI CI For 
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further details of this encoding, let us now have a look 
at examples. 



1. Encoding 2-body interactions, each site appearing in 
maximally n terms of the Hamiltonian function: Edge 
models 

A special case of the discussed generalization is the 
one where we (as in the preceding sections) stick to 
Hamiltonian functions with 2-body terms where each 
site is involved in n interactions. This kind of Hamilto- 
nian function plays a role in higher dimensional lattices 
and spin glasses, for example. For their treatment we 
propose, in the following, a way to discriminate the clas- 
sical spin configurations beyond resolving relative states 
(as we did before). 

To create a quantum state \jg) that enables us 
to differentiate energetically between all possible spin 
configurations, we proceed as follows. We identify 
two qudits with each edge e = (ij) G E of the 
graph G and provide them with a product basis 

{l s «) ei \ s j) ej \ s i> s j = 0...<7— lj, where is one of the 
edge qudits and ej is the other one. These qudits will 
be called the edge qudits corresponding to the edge. We 
map states of the classical spin sites to quantum state 
of the whole quantum many body system of edge qudits 
via 

Zf 38=(s„,..., S ,)h (g) \ Si ) ei \ Sj ) e _. 

e=(i])eE 

This way, we attach GHZ-states to the vertices, with 
the number of particles equaling the number of inci- 
dent edges. A graphical representation of this encoding 
is given in Fig. [3] (a). Note that each classical spin is 
mapped to as many edge qudits as there are edges at- 
tached to the classical spin vertex. 

Lemma 9. The superposition of the quantum states be- 
longing to all possible classical states 

l7G> == E ® l*>« l s i)e, 

seZ" e=(ij)EE 

is a product of GHZ-states and hence a stabilizer state. 

Proof. A reordering of the sites groups all edge qubits 
belonging to the same spin site i 



The overlap to evaluate the partition function or cor- 
relation functions has now to be performed with one 
state per edge qubit pair (a e \siSj). Since this overlap- 
ping state | a) allows us to adapt the energies hy to each 
individual spin, the possibility of evaluating partition 
functions with local energy terms as well as correlation 
functions is immediately given. 

To avoid the necessity of encoding the correlation 
function directly into the states \a> e ), we add one more 
quantum site to the GHZ state. This enables us to mea- 
sure the state of the classical site directly. Keep in mind 
that this is technically not necessary, because the state 
of the site is directly accessible already without the ex- 
tension. 

Theorem 10. The partition function Zq ({h e ,b v } ,f3) 
of a classical spin system at inverse temperature j3, de- 
fined on the graph G = (V, E) by the Hamiltonian func- 
tion H {{si}) = J2(ij)eE ( Si ' s j) ' can be written as 
the overlap of a stabilizer state and a product state ( over 
edge qudit 



wht 



Z G ({hi Jt b v },0) = ( <S> («(«)|)l7G>, 



Proof. The state \jg) is a product of GHZ states and 
hence a stabilizer state according to lemma and we 
compute, with an arbitrarily chosen orientation a of the 
graph G, 

( ® ( a (d)l7c> 

(ij)EE 

= ( <«w)l) E 

(ij)eE se z« e=(ij)eE 

= n e- phii ^^i) 



E 



-/3H({ Si }) 



□ 



\ S *)e t \ S j) ej 



e={ij)£E 



i£V e=(ij) 



and writing it this way we see that the state \ja) has 
the structure 

ho) = E ® K reo - ing (8) E K . 

sGZ™ i€V e=(ij) ieV SiGZ, e=(ij) 



where E Sie z 3 ® e = (ij) l s «)e; i s a GHZ state. 



□ 



2. Encoding A-body interactions, each site appearing in 
maximally 2 terms of the Hamiltonian function: Vertex 
models 

An important class of models are the vertex models. 
These models also fit into our framework, as will be 
shown now. A prominent example of a vertex model 
stems from a 2D regular lattice where each classical site 
interacts with two groups of three neighboring parti- 
cles (each individually) (see Fig. [4j). Hence we have 
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edge of edge of 




FIG. 3: Alternative encoding schemes I: Edge models. (The GHZ scheme.) This figure shows an example of an encoding of a 
classical interaction pattern into a product of GHZ states. The classical interaction graph, a square lattice in this example, 
is given by the underlying thin grid, the vertices symbolizing classical spin sites and the edges symbolizing their interactions. 
Each edge holds a pair of edge qubits, as indicated by the dots. The edge-qubits that belong to the same classical spin site 
are connected by thick lines, indicating that they form a GHZ state. The circle with dashed circumference indicates one pair 
of qubits \si) \sj) contracted with one state \ctij) in the Hilbert space of the pair of edge qubits. 



the situation where 4-body interactions take place, each 
site appearing in maximally 2 terms of the Hamiltonian 
function. Consequently, we encode each classical spin 
site into a 2-body GHZ-state (Bell state), and entangle 
the quartets of sites, corresponding to the interactions, 
in the state \a). 

This setting yields a vertex model in two dimensions, 
where the projections (of the subsets of the GHZ states 
taking place at each vertex of the vertex model) are de- 
termined by the set of states \a). Similar models in 
higher dimensions can be obtained easily in an analo- 
gous fashion. 



C. Relations between the encoding schemes 

An interesting question is how the generalized models 
that were just described relate to the encoding scheme 
encompassing the states \(Pg) and \iPg)- We want to dis- 
cuss this now and furthermore give additional relations 
between the states \<pc) and \iPg), adding to what was 
presented in the preceding sections. 



1. Relations between \<pg) and \iI>h) 

There are instances where for different graphs G and 
H the cutspaces of \<pg) and \iPh) are closely related. 
Two examples will now be demonstrated and give us 
some more insight into the internals of the construction. 

The first way to look at the construction of the 
cutspace of \<£g) is to modify the graph G by chang- 
ing the mapping of G to the quantum spin sites. We 
remember that (in the case of the construction of \iPg}), 
the method was to map each edge to one quantum spin 
site. As an alternative, we derive now from G an new 
graph by placing on each edge one additional vertex. 
This new graph we call the decorated graph G, which 
possesses N = \V\ + \E\ vertices. The crucial point is 
now to identify the vertices in G with the qudits that 
we chose as a product basis in the definition of the state 
\<Pg)- The original vertices (that appear in G and in 
G) are called vertex qudits and the qudits that were 
added at the edges are called edge qudits. The inci- 
dence matrix of the decorated graph G is now (ll-B ") 
with \<pg) = J2sez N I (^-\B a ) T s) , where for each sum- 
mand |s) \ (B a ) T s) the state |s) is a state of the vertex 

qudits and |(i? <T ) T s) is a state of the edge qudits. We 
note that the original method is a restriction of the just 
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i. . quantum 
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&-k4-c4>-cf 



measurement 
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(Bell pair ' 
in this case) 



\ v_7 / 
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FIG. 4: Alternative encoding schemes II: Vertex models. This figure shows an example of an encoding of a classical interaction 
pattern into a vertex model, where each thick line represents one Bell pair. In this example, each classical spin site enters 
in two four-site classical interactions. Accordingly, four edge qubits (circle with dashed circumference) form the smallest 
subsystem of the Hilbert space used for the overlap states \a £ ). 



proposed mapping of vertices to the edge qudits. 

A second way of mapping the graph G to another one 
that can be used to construct the cutspace of \<pg) is the 
following. Let us add one vertex to the graph G that is 
connected to all other vertices. Let us call this vertex 
h and the new graph G + h. The incidence matrix of 
G + h is 



B(G + h) J 



( B a (G) T \ 
1 



1 



V 



1 / 



The vector of classical spin sites s has to be extended 
to include the site h, hence we obtain a new vector s' = 
(s,Sh)- The canonical way to construct \ipG+h) now is 



G+h/ 



s' 

J2\B a (G) T s)\s + (s h ,s h ,...s h )) 

s,s h 

2]T|ir(G) T s )|s), 



because for all values of Sh, the equation 

B° (Gf (s + (s h , s h , ...s h ))) = \B a (G) t s 



holds. Hence \ipc+h) = 2|^g>- 

A conclusive remark seems appropriate. As has been 
shown, the stabilizer of the states are derived from the 
incidence matrix of their interaction graph. In the case 
of IV'g), the span of the rows of B a forms the cutspace 
directly. Following the arguments in the sections above, 
the stabilizer of \ipa) is constructed analogously, but 
from the span of the rows of the matrix (liyil-B' 7 ) or 

B (G + h) T instead. Although obviously being related, 
the difference in the construction changes the quantum 
states qualitatively to a great extend. For instance, 
when constructing {(fa) we do not obtain the same state 
classes as in the examples ©. Instead of the state 

(Sj=o (in case of a tree graph), or the state 

Sj^o 1 (\jx)® N ^J (in case of a cycle) or the toric code 
state, we always obtain states that are locally equivalent 
to one-dimensional or two-dimensional cluster states, re- 
spectively. 

Finally, by means of measurements, we are able to ob- 
tain the state I^g) from the state \ fa}- By overlapping 

the vertex qudits of \<pg) with the state (Xw=o U) ) 

we immediately recover ^g). On the one hand, this for- 
mally has the meaning of projecting out the dimensions 
of the state that are stabilized by operators correspond- 
ing to the l|y|-part in the matrix (l|y||B CT ). On the 
other hand, it has the physical interpretation of setting 
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the local external fields to zero. 



2. Going from the general picture to \<pg) and \iPg) 

The states \<pg) and ji/'c?) encode two-body interac- 
tions. Hence the scheme that is a direct super-set of the 
\<Pg) and \ipa) encodings is the GHZ scheme. Let this 
GHZ state be \GHZ). 

Contractions with states CK(#)) that do not discrimi- 
nate between quantum states |s;) e . \sj) e . with the same 
value of | Si — Sj\ yield directly the appropriate quan- 
tum description for \4>g}- To obtain the stabilizer de- 
scription for this case, each sub-state of \GHZ) con- 
sisting of a pair \si) e , \sj) e that is measured against a 

two-qudit state a^y)) is identified with a new single 
qudit. This qudit in turn corresponds to an edge in 
the adjacency matrix of the graph G defining the state 
\4>g)' The vertices G correspond to the GHZ sub-states 
in the state \GHZ). Hence all the information about 
the graph G can be recovered from the graphical scheme 
corresponding to \GHZ). The state \a!) that encodes 
the interaction strengths is not difficult to find either. 
Since Qf(y)) does not discriminate between quantum 
states \si) e . \sj) e . with the same value of |s< — Sj\ , we 
obtain 

K)} = E( E e-^-)W 

s \|a<-«jl 4 =» / 

Recovering a description of \GHZ) in terms of a 
state \<Pg) can be performed similarly, provided that 
the Hamiltonian function terms can be written as 

h(si, Sj) = hij (\si - SjQ + hi (s^ + hj (sj) . 

The GHZ state in the general encoding is a product 
state of smaller GHZ states 

\GHZ) = (g) \GHZ k ) . 
k 

Each of the states \GHZk) has to be extended by one 
site by the mapping 

N k N k + 1 

\GHZ k ) =^0| S ), ~ £ ® \s) t =: \GHZ' k ) . 

s i—1 s i—1 

To obtain the stabilizer description for this case, each 
sub-state of \GHZ) consisting of a pair |sj) e . \ s j) ej that 

is contracted with a two-qudit state is identified 

with a new single (edge) qudit of the decorated graph 
corresponding to \<pg)> This edge-qudit in turn corre- 
sponds to an edge in the adjacency matrix B of the 
graph G = (1\B) defining the state \<pa)- The sub- 
states that are not measured this way are the ones that 
were added in the mapping above. These will be used to 



encode the local fields and hence will be mapped to the 
vertex qudits of the decorated graph defining the state 
\<Pg)- The part that is more complicated here than in 
the case of \<£g) is finding the new state \a') encoding 
the interaction strengths. To do so, we have to find, for 
each term of the Hamiltonian function h (s,, Sj), a cor- 
responding form h (si, Sj) = hij (\si - Sj-M + hi (st) + 

hj (sj)- The part h (\si — s j\qj wm be encoded in the 
part of | a') that is measured against the edge-qudits, 
e.g., 

4i)) = E( E e-<Ml-«l.)V,). 

s \|si-sjl 9 =s / 

The local field corresponding to the vertex qudit with 
states \sk) Nk+1 , belonging to the extended GHZ sub- 
state \GHZk), is found by a summation of all corre- 
sponding fields 

N k 

h Nk+1 (s k ) = ^2hj(s k ), 
j=l 

where hj (sj) are the new terms of the Hamiltonian func- 
tion gained from the original terms h(si,Sj), that be- 
long to measurements on sites on the GHZ sub-state 
\GHZ k ). 

VI. APPLICATIONS 

This section contains applications of the framework 
given in the sections above. The first application shows 
how to derive the relation between a classical spin model 
on a graph and the corresponding model on the dual 
graph. The second application shows the implications 
of quantum mechanical symmetries existing in our de- 
scription of classical systems by means of a quantum 
system. 

Finally, we investigate the possibility to use simula- 
tions of the quantum system on a classical computer in 
order to obtain the statistics of quantum measurement 
results. This investigation yields some insight into the 
complexity of the computation of the partition func- 
tion and correlation functions of the classical system. 
We give a sufficient criterion for the structure of the 
interaction graph of the classical model, such that the 
computation of the partition function and correlation 
functions scale polynomially with system size. 

A. Duality relations for planar graphs 

We review Ref. [ll|. From graph theory it is known 
that for any planar graph G we can construct its dual 
graph D. In this section we want to demonstrate that 
the partition function Z G of a classical spin model de- 
fined on the graph G and the partition function Z D of 
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the model derived on the corresponding dual graph D 
have a simple and meaningful relation. 

To show this, we note that any orientation a of a 
graph G a induces an orientation of its dual graph D [lfjj , 
which we also denote by a (we refer to (l6fl . page 168 
for details). Moreover the incidence matrices B(D CT ) 
and B {G") corresponding to the two graphs have the 
property B (G a ) B (D a ) T = and the spaces gener- 
ated by the rows of these matrices are each others duals 
C-G (?) = Cjj (q) . Hence, the stabilizer of \iPd) can be 
written as 

S|. 0D) = {X (v) Z (u) \v eC D (q),ue C G (q)} . 
The quantum Fourier transform, 

9-1 



i * : 

v y jt k=o 



2-Kikj 

e q \j)(k\, 



has the property to map X and Z to each other un- 
der conjugation: FXF^ = Z and FZF^ = X, and can 
accordingly be used to map to <S|,/, G ), one-to-one, 

since F® N X (v) Z (u) (F® N )^ = Z (v) X (u) . Consider- 
ing the identity 



PS = 



g es 



for the density matrix ps of a stabilizer state that is 
stabilized by the q N operators in S, we infer that 

\1> D ) = F® N \ij G ) . 

The corresponding partition function Zq can thus be 
rewritten as 




where \a' e ) = F' |a e ). This transformation carries 
over to the energy terms in the Hamiltonian func- 
tion of the model on the dual graph, where we find 
Zg {q, <t, {h e }) = Zd (q, <r, {h' e }) with new energy terms 
h' e , which are derived from the old ones by 



-0KU) 



— Y< 



-0h e (k) 



for every j = 0, q — 1. 

We now want to examine the relation of the the Potts 
model on a graph G without external field and its corre- 
sponding model on the dual graph D. The Potts model, 
characterized by the Hamiltonian function 



H({ Si })=- J2 

e=(»,i) 



is encoded in two quantum states, \tpG) and &) e eE \ a e 
with 



a,, 



a 0J, 



9-1 

\o) + T\j) 



3=1 



The application of F^ on \a e ) yields 



q i/2 e -f3h' e (j) = I ' 
I e 



0Je + q -l if j = 
" J '-l if J = !,...,<?-!. 



Since the energies are again the same for all j = 1, q— 
1, we have another Potts model (on the dual graph D) 
whose interaction strength J' e fulfills the relation 



cfiJ'. - 



e PJ. +q -l 



Equivalently, we write (e 13 ^ — 1^ {eP Jc — l) = q, and 

hence recover the well known high-low temperature du- 
ality relation for the Potts model partition function [l|. 

B. Local symmetries 

See Ref. jljj . Local symmetries of stabilizer states 
can be used to show that several different models of 
classical spin systems actually have the same partition 
functions. More precisely, any local unitary U = ® e U e 
operator with eigenstate \ip G ) 



U |Vg) = A |Vg) 



(11) 



generates a model with the same interaction pattern 
but modified interaction strengths. Using Eq. l[TTj) we 
obtain the symmetry relation 

((g (a e |)h/> G ) = ((g) (OeDhfe) 



eeE 



where 



eeE 



The mapping 



9-1 

E 

3=0 



eeE 



\*U($\a e ). 

eeE 



9-1 



U e \a e ) =J2^ He(j) Ue\j) 
3=0 



implies another mapping of the energies defining the 
prefactors of the basis states This can lead to un- 
physical interaction strengths, e.g., imaginary ones. 

Similarly, a relation for the states \<pg) can be found, 
where the local symmetry is now corresponding to a 
change of interaction strengths and local field strengths 

((g K| (g (Oel) Wg) = ((g K|(g) (Oe|) \<p G ) , 

vev eeE vev eeE 
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where 

(g) K) (g) |5e) - A*t/ (g (g |Oe) . 

tieV e£_E v£V eeE 

The effect on the correlation function is again similar, 
but generically different correlation functions will, by 
the same symmetry transformation, be mapped to the 
corresponding correlation functions of different models. 
By definition, the state |a) enabling us to read out the 
value {s lll s l21 ...,s ln ) l3 is (g>„ ey \a' v ) ® e£is \ot e ) with 

|o.) = Ee-^W|i) 
j=o 

9-1 

\a' v (h, ...,*„)) = £ cos (27^)™" e-^W |j) , 
i=o 

where m„ is the number of occurrences of v in the n- 
tuple (ii,...,i n ). Now 

9-1 

L/„ |< *n)) = ^ £ cos (27ri/g) w " e"^ |j) , 

j=o 

so in general not only h e (j) and b v (j) will be al- 
tered, but the prefactors cos (2irj/q) m " play the role of 
weights. These are specific for the correlation function 
in question and enter the calculation of the energy terms 
b v (j) belonging to the symmetry. 

The fact that the states \i/jg) and \ipc) are stabi- 
lizer states is advantageous, because all elements from 
the stabilizer define such a symmetry operation already, 
which we will use in the following examples. 

Example 1 1 . We consider now the change of a classical 
model with q — 2 encoded into a state \ipc), caused by 
a symmetry operation. Let the classical graph have a 
vertex a with a set of edges E a connecting to it. One 
column c a of the incidence matrix corresponds to the 
vertex a. The stabilizer element X (c a ) Z (0) applied to 
the state \a) (encoding the interaction strengths) maps 
all interactions strengths J e ,e S E a to — J e and does 
not touch the other ones. We hence obtain the result 
that 

Z({J e }) = z({j e }), 

where 




- J e e e E a 
J e otherwise 



Next, we consider the change of a classical model 
with q = 2 encoded into a state \<pc), which is caused 
by a symmetry operation. The matrix generating the 
cutspace is now C — (1\B) . The construction of the 
local unitary symmetry operation using one column of 
C, like in the example above, yields now 

Z{{b u , J e }) = z({b„, J e }) , 



where 




- J e e S E a 
J e otherwise 



and b a = —b a and b v = b u otherwise. 



C. Simulations on classical computers 

An interesting aspect of the proposed mapping from 
classical to quantum systems is the established link 
between two different mathematical formalisms. As 
shown, algorithms for the computation of overlaps of 
stabilizer states with product states can be used to com- 
pute partition sums and correlation functions of classical 
spin systems - and vice versa. In both cases, hard and 
computationally feasible instances of these calculations 
are known, and we can now extend efficient algorithms 
from one domain to the other. This connection allows 
us to prove the following 

Theorem 12. There exists an algorithm that allows 
one to compute the partition function and the correla- 
tion functions of classical spin models defined on graphs 
exactly and with an effort that scales polynomially in the 
number of spin sites, provided that the tree-width of the 
graph used to define the classical model scales logarith- 
mically in the number of spin sites. 

The proof is rather technical and is given in ap- 
pendix El Thus, one finds that partition functions on 
graphs which are sufficiently similar to a tree graph (a 
property made precise by the notion of tree-width) can 
be efficiently evaluated. Similar results have been ob- 
tained in, e.g., Refs. [3. 



D. Relations to measurement based quantum 
computation 

In this section we discuss how the mappings between 
classical spin systems and the quantum stabilizer for- 
malism presented in this work, may provide insights in 
the study of measurement-based (or "one-way") quan- 
tum computation (MQC). 

The one-way quantum computer is a model of quan- 
tum computation introduced in Ref. [10]. In contrast to 
the quantum circuit model, where quantum computa- 
tions proceed by unitary evolutions, in MQC any com- 
putation is realized via single-qubit measurements only. 
More precisely, a one-way quantum computation essen- 
tially consists of two main steps: first, a system of many 
qubits is prepared in a highly entangled state, the "2D 
cluster state" [19], which is an instance of a stabilizer 
state. Second, part (possibly all) of the qubits in the 
system are measured individually. The qubits are mea- 
sured one after the other in a specific order, and each 
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qubit is measured in a certain basis which may (and typ- 
ically does) depend on the outcomes of previous mea- 
surements. It is this "measurement pattern" which de- 
termines the quantum algorithm which is implemented. 

It was shown in Refs. [jUllO] that the one-way quan- 
tum computer is a universal model for quantum com- 
puter, i.e., it is capable of (efficiently) simulating every 
quantum computation performed within the quantum 
circuit model. We refer to Ref. [13] for more details 
about MQC. 

Note that the model of MQC exhibits a remarkable 
feature, namely that the entire resource of a quantum 
computation is carried by the entangled cluster state in 
which the system is initially prepared. Indeed, as local 
measurements can only destroy entanglement, all the 
entanglement present within a one-way quantum com- 
putation must be provided by the initial resource state. 
Therefore, in order to understand the computational 
power of quantum computers, a study of the proper- 
ties of 2D cluster states, and other resource states, is 
called for. 

Even though it is by now well-established that the 
2D cluster states are universal resource states for MQC 
(and several other states have also been found to be 
universal [U, H2], it is not yet fully understood which 
properties of these states are responsible for their uni- 
versality. This issue has been the topic of recent in- 
vestigations [U, H3| (see also 0, E3, 0, El), where it 
was studied under which conditions a given quantum 
state may be a universal resource for MQC, and under 
which conditions it does not provide any computational 
speed-up with respect to classical computation. While 
significant progress has been made in these works, this 
important problem is far from being fully understood. 

What can the present connections between classical 
spin systems and quantum stabilizer states teach us 
about MQC? To this end, consider a one-way compu- 
tation having one of the stabilizer states \(pc) or \tpc) 
as a resource, where G is some graph. One may then 
ask which computational power can such resource states 
provide for MQC - i.e., which states among the \<pg) an d 
\tpa) are universal resource states, and which states are 
fully simulatable classically. Next we will see how the 
relation between these quantum states and the associ- 
ated classical spin systems, as established in this paper, 
provides insights in this issue. 

To do so, consider Eq. (!]), which identifies overlaps 
between a resource state \t)g) (= \ipc) or |</?g)) and a 
product state \a), as the partition function Zq of the as- 
sociated classical spin model on the graph G. Now note 
that such overlaps (to be precise, their squared modu- 
lus) equal the probabilities of outcomes of local measure- 
ments performed on the resource state \t]g). Therefore, 
if it is possible to compute such overlaps (and thus the 
corresponding measurement probabilities) efficiently, it 
becomes possible to simulate local measurement pro- 
cesses on such a resource, on a classical computer. Re- 
sources for which such efficient classical simulation is 
possible, by definition cannot offer any computational 



speed-up as compared to classical computation. Us- 
ing Eq. (I} , we now see that the problem of computing 
measurement probabilities of local measurements boils 
down to the evaluation of the partition function of the 
associated classical model. In particular, we find that 
classical models which are "solvable" — i.e., their parti- 
tion function can be efficiently evaluated — give rise to 
resource states for which the associated probabilities of 
local measurements can be computed efficiently. There- 
fore, the present mappings establish a relation between 
the solvability of a classical spin systems and the com- 
putational power of the associated resource state. 

Let us illustrate these relations with some examples 
for Ising models on different lattice types, with or with- 
out magnetic fields (see also Figs. Q] and [2]). Consider 
e.g., the simple case of a ID Ising model with periodic 
boundary conditions, without external field. This model 
is known to be solvable: its partition function can be 
evaluated in a time which scales polynomially with the 
number of spins. Using our correspondence, the associ- 
ated quantum state \iPg) is a GHZ state (see example 
[5]). This state is known to be an efficiently classically 
simulatable resource state for MQC. A similar conclu- 
sion can be drawn for the ID Ising model in the presence 
of an external field, which is solvable as well. Using our 
mappings, the associated quantum state \<Pg) is a ID 
cluster state, which is indeed also known to be simu- 
latable (see, e.g., 0]). Finally, also the 2D Ising model 
without field is known to be solvable - this is Onsager's 
famous result. The corresponding stabilizer state \4>g) 
is the toric code state. And indeed, this state is a sim- 
ulatable resource - in fact, the latter property has been 
shown in Ref. |13] by using the relation between this 
state and the solvable 2D Ising model. 

An Ising model which is not solvable is the 2D Ising 
model in the presence of an external field. In fact, the 
evaluation of its partition function is an NP-hard prob- 
lem. The corresponding stabilizer state is the 2D (dec- 
orated) cluster state. Interestingly, this state is a uni- 
versal resource for MQC. Therefore, we find that also 
in this case the computational difficulty of a classical 
model is reflected in the quantum computational power 
of the associated quantum state. 



VII. SUMMARY AND CONCLUSION 

In this work, we have displayed several mappings from 
Hamiltonian functions of classical spin systems to states 
of quantum spin systems. We map the interaction pat- 
tern given by the Hamiltonian function of the classical 
system to quantum stabilizer states and the interaction 
strengths as well as local field strengths to quantum 
product states. The overlap of these states yields the 
macroscopic quantities of the thermal states of the clas- 
sical spin system: the partition function and correlation 
functions at freely selectable temperatures (which are 
also encoded into the product states). 

The described mappings circumfere different classes 
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of admissible Hamiltonian functions. From the original 
and exemplary approach [ll[ suited for two-body inter- 
actions without local fields, we derive a more generalized 
mapping capable to yield correlation functions as well as 
to include local fields. Finally, we introduce a version 
capable to treat arbitrary Hamiltonian functions with 
n-body terms. Each of these mappings is interesting in 
its own right and offers an individual viewpoint and in- 
dividual aspects in the formal approach. The relations 
between the different mappings were investigated. 

We moreover gave several applications of the pro- 
posed mappings, namely: a simple derivation of the du- 
ality relation of a graph and its dual; a simple derivation 
of the impact of local symmetries of the stabilizer state 
on the classical model described by it; a constructive 
proof of a sufficient criterion for the possibility to effi- 
ciently evaluate of the thermal quantities of a classical 
spin system on a classical computer; and we discussed 
the relation of the computational accessibility of a classi- 
cal spin system with the power of a quantum computer. 
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APPENDIX A: A PROOF FOR THE GIVEN 
NUMBER OF STABILIZER ELEMENTS 



and a, b € Z^ 

(f(Xa + fib)= a{ci\\a + fib) = 

CieC B (q) 

A ^2 Ci{a\a)+fji ^2 Ci(ci\b) = Xcp (a)+fi(p (b) , 

c t eC G (q) CiEC G (q) 

by the linearity of the scalar product. The kernel 
of tp, ker(yj), is the set C G (q) ± because being a 
(orthonormal) basis {ci} is independent. The range 
of tp, ran ((/?), is the set C G (q), because for every 
w € C G (q) we have w = J2i an d V ( w ) — 
T^aecoiq) J2 3 ^j c i( c i\ c j) = w - The homomorphism (p., 
as any module-homomorphism induces an isomorphism 

1q /ker (<p) -^ran (<p) , 

which provides us with the formula 

\C G {q) | | ker Ml 

relating the number of elements in these sets. This im- 
plies 

q N = \C G (q) X \\C G (q)\. (Al) 

The number on the r.h.s. equals the number of the con- 
structed operators X (v) Z (u), which are, as a set, iso- 
morphic to 

{(c,s)\ceC G (q),seC G (q) X }. 
This concludes the proof. □ 



In section IIV Al we constructed the set of opera- 
tors X (v) Z (u) [see Eq. JH])] with v S C G (q) and 
u 6 C G (q) , where by construction C G (q) is the Z 9 - 
sub-module of that is generated by the rows of the 
incidence matrix B a . For this set to be a stabilizer of 
the single state \tp G ) it is necessary that \ip G ) is a fixed 
point of these operators (as already shown in the indi- 
cated section) and that is has cardinality q N . The latter 
point we show now. 

Lemma 13. The number of independent operators gen- 
erated by X (v) Z (it) [see Eq. JH])/ with v € C G {q) and 
u e C G (q) 1 - is q N . 

Proof. We note that the module and hence also all 
its sub-modules are free modules. Accordingly we can 
chose a basis, from which the modules or sub-modules 
are generated respectively. With the scalar product 
(•| • •) we construct an orthonormal basis {a} and with 
it the following mapping 

(p : —> Z^, w h-> ^ Ci{ci\w). 
This is a module-homomorphism, since for A, [i € Z g 



APPENDIX B: TENSOR TREE NETWORKS 
AND TENSOR TREE STATES 

We follow an approach of Shi, Duan and Vidal and 
consider the description of states in terms of a tensor 
network with tree structure 0, EH. We now want to 
give a short overview of fundamental definitions and 
theorems concerning these tensor tree states (TTS). 

1. Basic definitions 

The building block of a tensor network are complex 
d\ x d 2 x ••■ x d n tensors with elements Ai 1 i 2 ...i n . The 
number n is called the rank of the tensor A and the num- 
ber dk is called the rank of the index The maximal 
number d that the indices can assume, d = maxt dk, is 
called the dimension of the tensor. A summation over 
two indices ii and jv of common rank of two tensors 
and AM, 

A [r ' s] . 

= V A [r] A [s] 

Z_^, hi 2 — (ii=k)...i n 3x02 — OV=fe)— 3„i ' 
k 
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is called a contraction of the indices U and jy , A set 
of tensors together with pairs of indices that are to be 
contracted is called a tensor network. The maximal di- 
mension D of all tensors, D = max,4 d[A], is called the 
dimension of the network. Tensor networks can be rep- 
resented by graphs: the each vertex of the graph cor- 
responding to one tensor of the network and each edge 
corresponding to one pair of contracted indices. The 
indices to be contracted are referred to as internal in- 
dices and the other ones as open. The notation of graph 
theory carry over to the tensor networks, e.g., we talk 
about "subcubic" tensor trees. A tree graph (network) is 
called subcubic if each vertex (tensor) has degree (rank) 
1 or 3. The vertices with rank 1 are called leaves. 

It is possible to write the coefficients ^4 S of a generic 
pure iV-qudit state \ip) — J2 s As\s), where {|s)} is a 
product basis, as a contraction of a fixed set of ten- 
sors. Trivially, one tensor of rank N and a dimension 
equal to the number of states of the qudits is sufficient. 
In fact, representations for any graph-structure can be 
found, provided the rank of the internal indices being 
sufficiently large. Depending on the internal structure 
of the state to be represented, even representations with 
internal indices of comparatively small rank might be 
found, hence reducing the number of complex parame- 
ters representing the network. This displays the princi- 
ple that the more structure there is in the state, the less 
information is (potentially) needed to settle the remain- 
ing degrees of freedom. Conversely, any tensor network 
with N open indices can be used to define a pure N- 
qubit state. 

As an illustrative example, the tree depicted in 
Fig. [5] a) corresponds to the state 

M = A° A 1 A 2 A 3 k\ 

I / / j / j !l!2l3 tlSOSl '2S2S3 »3S4S5 I / ' 



Another well known example of subcubic tensor tree 
states are the matrix product states (MPS) with open 
boundary conditions. They have the simple form 

l MPS > =EE4o Sl 1*0*0 Kns 2 \s*) A\ i2S3 | S3 > 



Lemma 14. Calculating the overlap of a complete prod- 
uct state of N qudits with a subcubic TTS of dimension 
D over N qudits has a complexity of at most O (A^D 3 ) . 

Proof. Let the product state be 



and the TTS be \t). The calculation of (a\r) is a con- 
traction of a subcubic tensor network where the leaves 
are tensors with values (ai\l). In a subcubic tree, there 
is at least one tensor with at least two leaves attached. 
A contraction of this tensor with its attached leaves re- 
quires an effort of order D 3 . If this tensor has three 
leaves attached we are done. If not, this tensor will now 
be a leaf tensor attached to one other tensor and the 
tree will still be subcubic. As before, the tree will now 
have at least one tensor which has at least two leaves 
attached. We continue this procedure and because there 
are N — 2 tensors in the tree, we end up with an effort 
of the order ND 3 . □ 
Hence the contraction of a family of states over N 
qudits, whose TTN-description has a dimension scaling 
polynomially in N, with product states of the appropri- 
ate Hilbert spaces scales polynomially in N. 



3. Entanglement in TTS 

We want to state one more important result concern- 
ing the entanglement content of TTS with dimension 
D. Since the rank of the index corresponding to an 
edge connecting two tensors is limited by this number 
D, only D linearly independent combinations of states 
corresponding to the sub-trees attached to this edge are 
possible. Hence we have 

Lemma 15. The number of Schmidt coefficients of a 
TTS with dimension D in a bipartition of the qudits that 
corresponds to cutting exactly one edge in the graph of 
the corresponding TTN is limited by D. This Schmidt 
number can be reached. □ 



■ A i N sl-,s N \sn-is n ) 



APPENDIX C: A PROOF OF THEOREM [JJ 



2. Efficient scaling 

A subcubic TTN with N open indices (representing 
a subcubic TTS of N qudits) and dimension D depends 
on at most O (A\D 3 ) complex parameters. Thus, a fam- 
ily of states over N qudits whose TTN-description has 
a dimension scaling polynomially in A" allows for a de- 
scription with a number of parameters scaling polyno- 
mially in N. Concerning the contractions of TTS with 
product states we obtain the following result, (see also 
Refs. 0,113) 



We will prove theorem [12] using the machinery devel- 
oped in the sections IIVI and El The underlying idea of 
the proof is to map the classical spin problem (of finding 
the partition function) to the corresponding quantum 
problem (of finding an overlap), which is then solved 
by a simulation on a classical computer. To treat the 
simulation aspect, we need some results from the the- 
ory of tree tensor networks 0, 0] . We will use the tree 
tensor networks to encode the stabilizer states which 
are the images of the interaction patterns of the classi- 
cal spin systems. The necessary notation and theorems 
have been summarized in the appendix [B] With the 
language developed there, we reformulate theorem [121 
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|s > \s,) |s 2 > |s 3 > |s 4 > |s 5 > 

FIG. 5: A subcubic tensor network. The vertices correspond to tensors (circles) or physical sites (squares, "leaves") respectively. 
Edges indicate contractions over common indices. The bipartition of a state corresponding to the cut of a single vertex cannot 
have any Schmidt-rank higher than the rank of the connecting indices. The class of states generated by tensor networks 
covers all possible pure states provided that the dimension of the network is sufficiently large. 



1. Simulation complexity for the states \tpa) and 

\i>a) 

Theorem 16. For the states \ipc) and \<pc) (as defined 
above) a tree tensor network description can be com- 
puted with an effort growing polynomially in the number 
of classical spin sites N , provided that the tree-width of 
the graph G grows logarithmically in N . This tree ten- 
sor network description allows to compute the overlaps 
(a\ipG} and (a'\fG) of these states with product states 
with an effort that grows polynomially in N. 

Proof. The proof consists of three parts, i) In a 
preparatory step, we will summarize the ties between 
the tree-width of G and the branch-width of its cycle 
matroid. ii) We will then use this result to derive a 
bound for the Schmidt-rank of a TTN-description of the 
states \ipG) and |<Pg) and hence derive an upper bound 
for the computational effort to compute the overlaps 
{a\ipc) and (a'\ipG)- in) Finally we give an algorithm 
to find the (tensor-) coefficients in the TTN-description. 
Parts ii) and iii) have been given in a similar form for 
q = 2 already in Ref. Q. 

i) Let us first fix some notation, which can be found in 
more detail, together with missing proofs, for example 
in the references [H, EE El ■ 

Definition 17. (matroid) A matroid is a set f2 together 
with a rank function rk on its subsets. A rank function 
fulfills the following properties 

• If A and B are subsets of f2 and A C B, then 
rk(A) < rk(B). 

• For all subsets A and B of fi, 

rk (AnB)+ rk (A U B) < rk (A) + rk (B) . 

• If Ac fl, then rk (A) < |A| . 

A special instance of a matroid is the set of columns 
of the incidence matrix B a of a graph G — (V,E), called 
the cycle matroid M (G) of the graph G. 

A natural choice of a rank function on a cycle matroid 
is the dimension of the span of the column vectors. Fol- 
lowing the ideas of lemma [2] we deduce that with this 



choice of rank function and for a subset of column vec- 
tors A we have the relation rk (A) = \V\ — c, where c 
is the number of connected components in the graph 
G A = (V,A). 

Definition 18. (connectivity function) With the rank 
function rk of the cycle matroid M (G) of the graph 
G = (V, E) we define the connectivity function X on a 
subset of edges Ac E by 

A (A) := rk (A) +rk(E-A)- rk (E) + 1. 

It is a symmetric function with respect to A <-> E — A. 
An important observation is that with rk (A) = \V\ — 
c(A) follows the equality 

A (A) := \V\ + c(E) - c(A) - c(E - A) + 1, (CI) 

where c (E) , c (A) and c(E — A) are the numbers of 
the connected components in the respective subsets of 
edges. 

Definition 19. (branch decomposition) Let T be a sub- 
cubic tree (see appendix [B] and Fig. [5]) with edges E. 
The deletion of an edge e £ E of the tree corresponds 
to a bipartition of the set of leaves of the tree, because 
the deletion divides the tree into two (connected) com- 
ponents. The set of bipartitions of leaves induced by a 
tree is called a branch decomposition of the leaves. 

In the following we will identify the edges of our (dec- 
orated) interaction graph G describing the classical spin 
system with the leaves of a suitable tree. The set of pos- 
sible trees then corresponds to a set of different branch 
decompositions. With the connectivity function we can 
define the branch width of a branch decomposition. 

Definition 20. (branch-width) The branch-width 
br (A) associated with the branch decomposition in- 
duced by a tree T with edges Et is the value 

br (A) := max A (A e ) , 

where A e C E(T) is the subset of edges belonging to 
one of the remaining contiguous sub-trees of tree T ob- 
tained by deleting edge e £ E(T) from tree T [note the 
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symmetry X(A e ) = X(T — A e ).] The branch-width of the 
cycle matroid of G is defined as 



b (A) := min bx (A) . 



With this notation, we formulate the following theo- 
rems, to be found, together with the proofs, in Ref. (26| . 

Lemma 21. (Theorem 3.2 in Ref. fM]) Let G be a 
graph with at least one edge, and let M (G) be the cycle 
matroid ofG. Then the tree-width of G equals the tree- 
width of M(G). ' □ 

and 

Lemma 22. (Theorem 4.2 in Ref. fM]) Let M be a 
matroid of tree-width t and branch-width b. Then 



Lemma 24. Let PUQ = E with PnQ = be a biparti- 
tion of the edges in the interaction graph G = (V, E) de- 
scribing the interaction pattern of the classical spin sys- 
tem. Let \tig) denote the quantum state whose stabilizer 
is constructed via M (e.g. \(p G ) or \^ G )). The Schmidt 

rank \ °f the bipartition \rj G ) = Y^f=x I^D 



b-l<t< max(2b-2, 1) . 



□ 



In particular this result tells us that the tree-width 
is an upper bound to the branch-width. The next the- 
orem, to be found in Ref. [23], now states that we can 
algorithmically compute a subcubic tree that at least 
comes close to the optimal tree. 

Lemma 23. (Theorem 2.12 in Ref. \21i]) For given 
k, there is an algorithm as follows. It takes as input 
a finite set Eg with \E G \ > 2 [and the connectivity 
function X]. It either concludes that 6(A) > k or out- 
puts a tree with br (A) < 3k + 1. Its running time is 

O 

□ 



where | f]f ) and rff \ are quantum states of the qudits 

corresponding to the edges in P and Q respectively, sat- 
isfies the equality 

y = q \V\+c(E)-c{P)-c{Q) 

where c (E) , c (P) , c (Q) are the number of connected 
components in the graphs (V, E) , (V, P) and (V, Q) re- 
spectively. 

Proof. Corresponding to the bipartition P U Q = E 
we have a bipartition of the columns of the matrix 
M. After performing some (unimportant) permuta- 
tion of the columns, the matrix M takes the form 
M = {Mp\Mq). Let c denote the number of columns of 
M and p and q the number of columns of Mp and Mq 
respectively. 

Now let S be the stabilizer of the state \vg) and Sp c 
S the subset of operators g that act trivially on the 
qudits belonging to the labels in Q. We define Sp := 
{Tr<2 Sp}. From the theory of stabilizer states 

it is known that |t?g)(77g| = q~ c J2 g es9 an( ^ nence 



Pp 



^|S G | 6 log|£; G ^, where 5 is the time to compute X. Tr Q [\7] G )(n G \] = Ttq 



Because efficient algorithms to compute the tree- 
width of a graph G (and hence with the lemma above, 
upper bounds for the branch-width) are known, we can 
assume to be able to input a k > b (A). This way we al- 
ways end up with a tree T with at most br (A) = 3fc+ 1. 

At the end of this part we know, given an interaction 
graph G with tree-width t, that we can efficiently com- 
pute a branch decomposition over the set of edges such 
that the branch-width associated with this decomposi- 
tion is smaller than or equal to 3t. 

ii) We now want to establish a link between the 
X~ width associated with bipartitions of the states \iPg) 
and \^pg) and the branch- width of the a tree inducing 
these bipartitions. To fix some notation, we define the 
matrix 



g&Sp 

= T P E 9- 

gESp 



The factor g q comes in because the trace over all oper- 
ators but 1 in the Pauli group is zero and Trg[l] = q q . 
Furthermore, the stabilizer is a group, so we have the 
identity 



(Pp) = 



E E^ 



s F 



,2p 



E 



heSf 
\Sp\ 



-pp. 




for \tp G ) 
for \ip G ) 



(C2) 



We recognize that M is used to define the stabilizer of 
the respective states, because 



C G = {Af T s,s e . 



We define r :— q v / \Sp\ and obtain (rpp) = rpp. 
Hence rpp is a projector and has (after a possibly 
necessary change of basis and reordering of rows) the 
form rpp = diag (1, 1, 0, 0), or equivalently, pp = 
diag (r _1 , r _1 , 0, 0). Since Tr[pp] = 1, we have 
r -1 rank(/5p) = 1 and hence r equals the number of 
Schmidt coefficients in the bipartition of the state \r) G ) 
according to the sets of edges P and Q. Thus X = r = 

q p /\s P \. 
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To obtain the number \Sp\, we have now a look at 
the matrix M — (Mp\Mq), which we will from now on 



interpret as a linear mapping M 1 : Zg V ' — > IS Here, 
Mp is a |V| x p-matrix belonging to the columns in P 
and Mq is a \ V\ x q-matrix belonging to the columns in 
Q. Recall that the stabilizer is isomorphic to the set of 
operators 



X(v)Z(u) 



of M 



where v E Cg (q) and u E Cg (q) . Hence |<Sp| is de- 
termined by the number of vectors v' C Cq (q) and 
u' € Cq (q) whose elements are in the last q places 
(e.g. v' = (v[, v' p , 0, 0)). Let this number for the 
set C G (q) be z c = \C P (q)\, where (v[, ...,v' p ) E C P (q), 

and the corresponding number for the set Cg (q) ± be 
z c ± = \Cp-(q)\, where (u[, ...,u' p ) S Cp-(q). Then 
\S P \ = z c z c ±. 

Let us now calculate zc ■ The elements of Cg are the 
image vectors of M T . Furthermore, if s S ker(-Mg), 
then the image of s has the desired form M T s = v' = 
(v[, —,vL, 0, 0). Considering that we can add any 
vector from the kernel of M T to s without changing the 
image v, it is zc = |ker {Mq) \ / |ker (M T ) | . 

Similarly, z c ± equals the number of elements in the 
set Cp(q) 1 ' where Cp{q) — ran(Mj). As shown as 
part of Appendix A, this number is equal to z C -l = 
q p I ran (Mp) | as the target space of the mapping Mp 
isZS 



Another basic consideration about the linear mapping 
v\ 

spaces) tells us that 



Ml 



7L}a 1 — > ran (Mj) (note: a mapping between finite 



hence 



and 



q 



zczc^ 



m = Iran (M£)| |ker(M?) 



\ker (Ml) 



\ker (M%)\ 



\ker (M T )\ 



q W\ 


ker (M T ) 




fcer(Mj) 




ker (Mq) 



From lemma [2] we now derive that \ker (M T ) 
(with analogous results for Mp and Mq). 

Remark 25. As a side remark we note the identities 



□ 



r = q»/\Cp{q)\\Cp-(q)\ 



(C3) 



and 



|ran(Mj)| = \C P (q)\ = q^ / |ker (Mj) | - q \ v \- c ( p \ 
which can be obtained from the proof above. 



Corollary 26. Considering identity fCij) . we deduce 
that the Schmidt-rank \ of a bipartition of edge-qudits 
E = AlJ (E — A) and the connectivity function A defined 
on the graph G satisfy the following equation 



X 



,A(A)-i 



Considering that the matrix M defined in Eq. I|C2|1 is 
the cycle matroid linked to the states (i/'g) and \(pc) we 
can now state the following important result, concluding 
the second part of the proof. 

Corollary 27. Using the result of lemma HH to find, 
by means of the matrix M , a branch decomposition of 
the audits in the states \iPg) and \ f>G), uie can efficiently 
find a subcubic TTN description such that the Schmidt 
number of all bipartitions following this branch decom- 
position satisfies 



X. 



< 



According to lemma \TE[ the dimension D of this TTS 
is limited by it — 1 and hence, following lemma\Tl\ the 
effort to compute the overlaps (a\ipG) and (a'\(pG) grows 
with at most O (\E G \t 3 ). □ 

iii) In this part we want to discuss how to compute 
the tensor entries in the TTS description of the states 
\ipc) and \<pa), which we will again denote generically 
as \r)c) where no distinction is necessary. The ansatz 
for the calculation of all tensor elements is the branch 
decomposition of the edge qudits (concerning the edges 
of the graph G = (V,Eg)) induced by the tree tensor 
network T with edges Et describing the state. We select 
an arbitrary edge eg € Et of the tree to obtain an initial 
bipartition Eg = P^Q with Q = (Eg — P) of the edges 
in Eg , inducing a bipartition of the set of qudits of the 
state \t)g) ■ We will use the notation P and Q for the 
edges and the corresponding qudits alike. 

Let us consider the Schmidt decomposition belong- 
ing to the bipartition. Recalling the proof of lemma 
l24l the Schmidt coefficients of a decomposition \t)g) = 
Tli \Pi) \qi)i where the states \pi) live on the Hilbert 
space of the edge qudits in a part P C Eg and the 
states \qi) live on the part Q — Eg — P C Eg can be 
obtained immediate^ They are all equal and have the 
value Aj = r _1 = |<Sp| /q p . We remember also that 
there are exactly r of these coefficients. Concerning the 
Schmidt vectors, we consider the following lemma. 

Lemma 28. A Schmidt basis for a bipartition of the 
edge qudits Eg = P U Q, Q = Eg — P of the state \tjg) 
is given by the set of states {\pi) |<7z)}i=i where 

\ Pi) := q ip-\ym ^ \ Cp+Pi) 

cpGCp 

h) := ^- |V|)/2 £ \c Q + q i ). 

cqECq 

Here p~i E (Cp-) 1 " , such that the cosets Pi + Cp are dis- 
tinct for different values of i, and Cp is the cut space of 
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the subspace belonging to the edges belonging to the edge 
qudits in P. ({q~i} C (Cq) is defined analogously; all 
additions in the kets are modulo q). 

Proof. We look at the states \pi) first; the states 
\qi) are treated analogously. The set of states {\pi)} 
has to be an orthonormal set which at the same time 
is a set of eigenstates of the reduced density operator 
PP = Trg [\vg)(Vg\]- We define S P := {Tr Q [g],ge S P } 
and recall from the proof of lemma l24l that |t7g)(?7g| = 
1~ c 12ges9 and hence p P = q~ p J2 g es P 9- Now each 
g £ Sp can be written as g = X (v) Z (u) where v £ Cp 
and u £ Cp. Applying such an operator to \pi) yields 



By construction, the vectors qi are all elements of 
(Cq) . Furthermore we notice that there are efficient 
algorithms for all these steps. □ 

This bipartition enables us to compute all tensor en- 
tries efficiently. 

Consider that using the TTN description of the state 
\t)g) the states \pi) and \qi) can be written as 

\pi) = E 41 ii>« i fc )p 2 > \pi) = E A ul h Q2 
jk 



-(P-IVD/2, 



\ Pi ) = X(v)Z{u) \CP+Pi) 

cpECp 

c P eCp 

= E I c p+^)' 

c'ecp 



since pi £ (cp) and cp is a group. The perpendic- 
ularity property of the states \pi) stems from the fact 
that the vectors pi are from distinct cosets for different 
values of i. We furthermore calculate 



(Pi\Pi 



q 



(P-\V\) 



E 

cp,c'£Cp 



6 cp:C' p = 1 



following from remark [25l The number of Schmidt vec- 
tors is indeed r, because the number of distinct cosets is, 
with a slight generalization of the results of Appendix 

lAl especially Eq. (|Al[) . to Cp and 



{<%) 



Cp-t/Ci 



CI 



(I 



\p\ 



Eq. lAll 



\c F 



rk[25l 



Having proven that the individual states \pi) and \qi) 
have the given form, we note that the pairing (pi,qi) for 
each i is not arbitrary and has to be found out. In the 
following we give an algorithm to find these pairs. We 
assume in this context that joining the edges of P and 
Q results in a mere concatenation of the corresponding 
vectors to simplify the notation. This can always be 
achieved by a reordering of the edges. The algorithm 
that we use is as follows 

1. Find the set {pi}, an orthonormal basis {cq} of 
the space Cq and a vector (cp|0) £ C, where cp £ 
C P . 

2. For each pi find one vector (cp + pi\a,i) £ C where 
the choice of is in principle arbitrary and just 
limited by the set of vectors in C. Keep the vectors 



3. For each vector a,i calculate the corresponding vec- 
tor ^ := a, - J2 5q 5 Q ( 5 Q ' a i) ■ 



jk 



with suitable tensors v4i p l and A™ and states 
\j) P , \k)p 2 . The states \j) P are living on the Hilbert 
space Pi belonging to the leaves (and hence to the 
corresponding qudits) that are part of the sub-tree 
of T attached to the tensor A^ by its index i. 
Also the states \j) p can be written as 



[Pi] 



/Pi 



>Ptt 



>Pl2 



(analogous arguments apply to 



the states \k)p 2 , \1}q 1 , \ m )Q 2 -) To be able to compute 
the entries of the tensors we hence need the states be- 
longing to sub-trees which can be derived from the ini- 
tial Schmidt decomposition. 

Lemma 29. Let \i) E be a state on the qudits cor- 
responding to a set of edges E, defined as \i) E 

E 



\ce + d(i)), where Ce is the cut space of the 
incidence matrix of the graph G — (V, E) belonging to 
the qudits as defined above. Let E ~ PUQ, Q = E — P 
be a bipartition of the qudits. Then 



whe 



V) E = \Pi) \Qi) 



\Pi) = E \tP + d{i) t 

c P £Cp 

\Q*) = E \ c Q + d ( l )c 

cqGCq 



The states |Pj) and \Qi) live on the Hilbert spaces of P 
and Q respectively and the vectors d(i) p and d(i)g are 
the parts of the vector d (i) belonging to the respective 



Proof. A reordering of the position of the qudits in 
\i) E , so that the merging of the vectors cp,CQ,d(i) P 
and d(i)g becomes a concatenation yields 



E 



\cQ) + (d(i) P \d(i) c 



cpECp.cqGCq 

E \(Cp\CQ)+d(i)). 

cpeCp.cgeCQ 

The sets of vectors {ce},{cp} and {cq} are the im- 
ages of the matrices M E ,Mp and Mq of identity IC2l 
respectively, where the index denotes the edges that the 
columns correspond to. Since Me = (M p \Mq), we ob- 
tain the result that for each ce there is exactly one pair 
(cp, cq) immediately. □ 



25 



We observe that the involved sets Cp,Cq and the 
vectors d(i) p ,d(i)g can be found efficiently. Now we 
write \i) E = J2i 3 k Aijk \Pj) \Qk) and deduce that 

except for AM and AM which have to absorb the square 
root of the Schmidt coefficients also and hence A^'^ = 
r- 1 / 2 S lJ S tk . 

This concludes the proof of theorem HH □ 

2. The bipartite entanglement of the general 
encoding schemes (e.g. GHZ-product state and the 
vertex model state) 

So far we have only considered the computational 
complexity using an encoding into the states \i^g) and 
\ific)- In this section we want to extend the efficiency 
statement to the alternative encoding schemes discussed 
in section Ivl 

The major modification leading to these schemes and 
complicating the entanglement aspect is the extension 
of measurements from one qudit to two or more. In a 
branch decomposition, the sites being involved in these 
measurements have to be placed in their own sub-trees, 
which we will refer to as contraction sites. The con- 
traction of the highly entangled states \a e ) with these 
contraction sites will in general not be efficient, but since 
the size of the contraction sites is limited, this only leads 
to a constant computational overhead. In a branch de- 
composition of a state of the extended encoding schemes 
we can represent the contraction sites as leaves. 



The remaining question is "What is the entanglement 
of bi-partitions in a branch decomposition where the 
contraction sites are leaves"? By construction, we im- 
mediately find that this question can be answered by 
looking at the number of states (in our schemes those 
are either q-dimensional Bell pairs or GHZ states) that 
are shared by different contraction sites and cut by the 
branch decomposition. See also Figs. [3] and [H 

Once we contract the contraction sites in the graph- 
ical representation of the general picture (like given in 
Fig. [3]) to single vertices, we obtain a new graph where 
the edges represent Bell pairs shared by contraction 
sites. Graph theory immediately tells us that also in this 
case the tree width is the decisive quantity of the (con- 
tracted) graph that governs the minimum number of 
states (and hence ebits) that have to be cut in a branch 
decomposition. The tree width of the contracted graph 
is carried over from the underlying graph of the classical 
interaction graph. Thus we can conclude that theorem 
[16] applies for the alternative encoding schemes as well, 
and the decisive parameters can be derived immediately 
from the respective encoding patterns. 

We also emphasize that non-planar graphs of loga- 
rithmically bounded tree-width, as well as non-local in- 
teractions are covered by this result. Results regarding 
efficient computation of homogeneous Potts model par- 
tition functions on graphs of logarithmically bounded 
tree-width have been obtained before, though with en- 
tirely different methods. We emphasize that our ap- 
proach, in contrast to previous approaches, can handle 
without difficulty also inhomogeneous models. More- 
over, it leaves a lot of space for generalizations. 
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